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Abstract 

The electronic spectrum, energy gap and local magnetic moment of paramagnetic NiO are com- 
puted by using the local density approximation plus dynamical mean-field theory (LDA+DMFT). 
To this end the noninteracting Hamiltonian obtained within the local density approximation (LDA) 
is expressed in Wannier functions basis, with only the five anti-bonding bands with mainly Ni 'id 
character taken into account. Complementing it by local Coulomb interactions one arrives at a 
material-specific many-body Hamiltonian which is solved by DMFT together with quantum Monte- 
Carlo (QMC) simulations. The large insulating gap in NiO is found to be a result of the strong 
electronic correlations in the paramagnetic state. In the vicinity of the gap region, the shape of 
the electronic spectrum calculated in this way is in good agreement with the experimental x-ray- 
photoemission and bremsstrahlung-isochromat-spectroscopy results of Sawatzky and Allen. The 
value of the local magnetic moment computed in the paramagnetic phase (PM) agrees well with 
that measured in the antiferromagnetic (AFM) phase. Our results for the electronic spectrum and 
the local magnetic moment in the PM phase are in accordance with the experimental finding that 
AFM long-range order has no significant influence on the electronic structure of NiO. 



PACS numbers: 71.27.+a, 71.30.+h, 79.60.-1 



I. INTRODUCTION 



NiO is a strongly correlated electron material with a large insulating gap of 4.3 eV and an 
antiferromagnetic (AFM) ordering temperature of Ti^=523 K. Conventional band theories 
which stress the delocalized nature of electrons cannot explain this large gap and predict NiO 
to be metallic- On the other hand, spin-polarized band calculations, e.g., density functional 
calculations based on the local spin-density approximation (LSDA),- which do find an AFM 
insulating ground state in NiO, produce a band gap and local magnetic moment which are 
considerably smaller than the experimental values. These facts are often taken as evidence 
for the inapplicability of conventional band theories to strongly correlated systems like NiO. 
Indeed, already a long time ago Mott^ had shown that NiO and similar insulators may 
be better understood within a real-space picture of the solid, where localized electrons are 
bound to atoms with incompletely filled shells. This leads to the formation of incoherent 
bands, the lower and upper Hubbard bands, which are separated by a correlation gap of 
the order of the local Coulomb repulsion U. For this reason NiO has long been viewed as a 
prototype "Mott insulator" .-i^ 

This view of NiO was later replaced by that of a "charge transfer insulator" ,- after Fu- 
jimori and Minami had successfully explained the photoemission data in terms of a cluster 
model treated within the configuration-interaction method.- In particular, this interpre- 
tation was supported by the combined x-ray-photoemission (XPS) and Bremsstrahlung- 
isochromat-spectroscopy (BIS) measurements of Sawatzky and Allen.- Within this new pic- 
ture an additional ligand p band appears between the lower and upper Hubbard bands, and 
the insulating gap is formed between the ligand p band and the upper Hubbard rf-band. 
However, unless the p-d hybridization is taken into account, this picture is still an oversim- 
plification. Namely, the hybridization between transition-metal d and ligand p states will 
lead to some (i-electron features also in the upper valence bands. Indeed, subsequent stud- 
ies suggested that the first valence peak is actually a bound state arising from the strong 
hybridization of Ni 3d and O 2p states,-*^ such that NiO is close to the intermediate regime 
of the Zaanen-Sawatzky-Allen scheme.- Despite the success of the cluster approach it has 
apparent drawbacks since it neglects the band aspects of O 2p states completely which are 
known to play an important role in NiO.-ii^ The translational symmetry has been taken into 
account to some extent within the cluster perturbation theory recently.— Another extension 
is the treatment of a larger cluster (NigOig) so that nonlocal charge transfer excitations can 
be identified.— 

Since the cluster approach relies on adjustable parameters to fit the experimental spec- 
trum it is highly desirable to obtain a description of the electronic structure of NiO from 
first principles. Already within L(S)DA the O 2p bands can be accounted for quite well.- 
Attempts to go beyond L(S)DA are based on the self-interaction-corrected density func- 
tional theory (SIC-DFT),— the LDA+U method,— and the GW approximation.— These 
methods represent corrections of the single-particle Kohn-Sham potential in one way or an- 
other, and lead to substantial improvements over the L(S)DA results for the values of the 
energy gap and local moment. Within the SIC-DFT and LDA-I-U methods the occupied 
and unoccupied states are split by the Coulomb interaction U, whereas within LSDA this 
splitting is caused by the Stoner parameter /, which is typically one order of magnitude 
smaller than U. Therefore, compared with LSDA, the SIC-DFT and LDA+U methods cap- 
ture more correctly the physics of transition-metal (TM) oxides, and improve the results for 
the energy gap and local moment significantly. The GW method goes one step further by 
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calculating the self-energy to lowest order in the screened Coulomb interaction W, and the 
obtained band structure shows better agreement with angle-resolved photoemission spectra 
(ARPES).— Since then a large number of studies of NiO have been performed along these 
lines,— ii^ii^'^i^>22i>2^ differing in the basis used and/or the details of the approximations. 
However, in both SIC-DFT and LDA+U the self-energy is static, i.e., does not take corre- 
lation effects into account adequately and thus cannot give an accurate description of the 
electronic energy spectrum of this correlated material. As to the GW method, its practi- 
cal applications to NiO indicate that, in general, it also cannot explain strong correlation 
effects sufficiently. Furthermore, different implementations of the GW scheme lead to quite 
different result s^'i^'^Si'^ regarding the value of the insulating gap and the relative positions 
of the energy bands. 

Thus, although considerable progress was made in the theoretical understanding of NiO 
from first principles, several important issues are still open: first of all, the agreement 
with electronic energy spectra is far from satisfactory. This is not really surprising since 
the self-energies employed in the previous approaches are either energy-independent (i.e., 
genuine correlations are not included) or are obtained from approximations whose validity 
is not entirely clear. Secondly, since conventional first-principles calculations are not able to 
incorporate the strong local electronic correlations they have, with only a few exceptions,— 
focused on the AFM ground state, and hence attributed the insulating gap to the existence 
of long-range magnetic order. It is well known, however, that both the band gap and the 
local magnetic moment are essentially unchanged even above the Neel temperature.- Indeed, 
recent experiments showed that the long-range magnetic order has no significant influence 
on the valence band photoemission spectra^^ and the electron density distributions.— These 
facts are evidence for the strongly localized nature of the electronic states in NiO and 
cannot be understood within a Slater-type, static mean-field theory of antiferromagnetism.— 
Apparently the AFM long-range order itself is not the driving force behind the opening of 
the insulating gap but is merely a concomitant phenomenon, i.e., a consequence of the 
correlation-induced insulating gap, rather than its origin. We are thus led to the question 
whether the properties of the strongly correlated, paramagnetic (PM) insulating state of 
NiO can be calculated within a theoretical framework, i.e., explicitly explained, even in the 
absence of long-range AFM order. 

In this work, we therefore investigate the PM insulating phase of NiO using the 
LDA+DMFT approach. This computational scheme for the ab initio investigation of 
strongly correlated materials was developed during the last decade^^i^Ei^ and has led to 
important insights into the physics of strongly correlated electron materials; for reviews see 
Refs. 3l| S^ SSjli 3^. The LDA+DMFT approach combines band structure theory within 
the Local Density Approximation (LDA) with many-body theory as provided by dynamical 
mean- field theory (DMFT).— s2I Within DMFT, a lattice model is mapped onto an effective 
impurity problem embedded in a medium which has to be determined self-consistently,— 
e.g., by quantum Monte-Carlo (QMC) simulations.— This mapping becomes exact in the 
limit of infinite dimensions.— 

The LDA+DMFT approach is a particularly suitable scheme also for the investigation 
of the electronic properties of NiO. Namely, it takes into account both the material specific 
aspects as well as the strong electronic correlations in NiO, thus having advantages that the 
previous methods lack. This approach has been used by Savrasov and Kotliar— to study the 
phonon spectrum of NiO and MnO. In this work we concentrate on the electronic spectrum. 

The paper is organized as follows. In section [TTl the LDA+DMFT computational scheme 
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for NiO is presented. The general idea of modeling NiO by a material-specific many-body 
Hamiltonian is discussed in section III A[ In section III Bl the Wannier functions and single- 
particle Hamiltonian matrix are constructed within the LDA, and the correlations are in- 
cluded by DMFT(QMC) in section III CI The results are discussed and compared with 
experiment in section IIII| and finally we conclude this paper in section IIVI 

II. COMPUTATIONAL METHOD 

A. Description of NiO with a multi-band Hubbard model 

In practice, LDA calculations usually involve a large number of valence s, p, and d 
orbitals associated with all the atoms in the unit cell. Since multi-orbital QMC calculations 
are computationally expensive, especially at low temperatures, not all of these orbitals can 
be included in the DMFT calculation. For this reason one needs to project out most of the 
orbitals except for the most relevant ones. 

For a large number of transition-metal compounds, the most relevant orbitals responsible 
for the physical properties are the transition-metal valence d orbitals which are partially 
filled and are located around the Fermi level. However, for some materials, e.g., the late 
transition metal oxides, which are either charge-transfer insulators or in the intermediate 
regime of the Zaanen-Sawatzky-AUen scheme,- the oxygen 2p orbitals are as important as 
the transition-metal valence d orbitals, and hence have to be taken into account in the 
effective model Hamiltonian as well. Thus, in principle, the strong Coulomb interaction 
among the d electrons and the p-d hybridization should be included and investigated on the 
same level. In DMFT this would correspond to the investigation of a material-specific multi- 
orbital version (five d orbitals and three p orbitals) of the periodic Anderson model with 
nearest-neighbor hybridization. Such a comprehensive DMFT treatment is very difficult to 
perform even at present. Therefore, as a first step we turn to an approximate, but useful 
procedure: namely, we include the hybridization only implicitly in the construction of the 
set of Wannier functions as will be described below.— With this set of Wannier functions as a 
basis, we arrive at an effective multi-band Hubbard model. The on-site Coulomb interaction 
for this Hamiltonian is then an effective interaction Ues for Wannier orbitals which is different 
from the bare Coulomb interaction U between the d electrons. In this way we are actually 
mimicking the charge transfer gap by a Mott gap. Such a simplified treatment is quite 
analogous to the reduction of the two-band p-d hybridized model to a one-band Hubbard 
model in the context of cuprates.— (We note in passing that the Hubbard model and the 
periodic Anderson model with nearest neighbor hybridization display a striking similarity 
near the metal-insulator transition.—) Our simplified treatment of the oxygen 2p bands 
and the p-d hybridization should be regarded as a first step towards a full LDA+DMFT 
investigation of the NiO problem. Within this approximation some of the spectral weight is 
lost and, in particular, the satellite structure at higher binding energies is ignored. Thus we 
are not aiming at a calculation of the full spectrum, but rather concentrate on the spectrum 
near the gap. 

Wannier functions can be constructed in different ways. Historically there were 
attempts^i^'^ to produce Wannier functions by a variational procedure, without know- 
ing the actual Bloch functions. Within the context of LDA+DMFT, Andersen and Saha- 
Dasgupta4^ proposed a way of constructing Wannier functions a priori based on A^th-order 
muffin-tin orbitals (NMTO) which was then applied to 3d^ perovskites^ and Ti203.— On the 
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other hand, Wannier functions can be calculated straightforwardly by a unitary transforma- 
tion of a set of Bloch bands. However, due to the nonuniqueness of Wannier functions, quan- 
titative computations along this line did not appear until recently. An important progress 
was made by Marzari and Vanderbilt^ who devised an optimization procedure to obtain 
the maximally localized Wannier functions. However, the maximal localization of Wannier 
functions is not a strict requirement if one is not interested in the orbitals themselves but 
in a proper representation of the Hamiltonian, which is the case in this work. 

In our paper, we follow a procedure proposed recently by Anisimov et al.,— in which 
a set of d-like Wannier functions^ is constructed in such a way that they (i) preserve 
the symmetries of the atomic-like d orbitals, and (ii) implicitly incorporate the admixture 
of oxygen 2p orbitals arising from the hybridization effects. The LDA band structure is 
encoded in a non-interacting Hamiltonian within this basis set of Wannier functions and 
used as the input for the DMFT study. In the above-mentioned work,— calculations were 
carried out for SrVOs and V2O3. Here we extend this scheme to NiO. 

From the LDA+DMFT calculation, one can obtain the energy gap, the local magnetic 
moment, and the electronic spectrum. Since these quantities do not significantly depend on 
temperature we compare our results with experimental data at low temperatures and other 
theoretical results for the ground state. 



B. LDA band structure and construction of Wannier functions 

We first perform a standard LDA band calculation using the linear muffin-tin orbital 
(LMTO) method in atomic sphere approximation^ method with the combined correction 
term included. The Stuttgart LMTO code version 47 is used for this calculation. Para- 
magnetic NiO has rock salt crystal structure, and below T/v there is a tiny rhombohedral 
distortion of the cubic unit cell. The lattice constant shows a small increase from about 
4.17 A to 4.20 A in a temperature range from 7 K to 700 K.— Such an increase only causes 
a small deviation of the bands relatively far away from the Fermi energy, but no noticeable 
change of the Ni 3(i-dominant bands on which we will focus later. Since we will compare our 
result with the low-temperature experimental data, we choose the lattice constant a = 4.17 A 
throughout this work. The calculated nonmagnetic band structure of NiO is shown in Fig. [1] 
and it is in agreement with those published in the literature.-'^i 

Based on the solution of the LDA band problem, we start to construct a set of 
nonorthogonal Wannier functions centering at a lattice site T, following Ref. [s^ . 

N2 

= ^ E^*''!^"'^) = ^ E^^''^ E \^^>^)i^^''\^n)- (1) 
^■^^ k ^ k i=Ni 

Here \4>ik) are the single-particle Bloch states and are the LMTO basis functions. 
Furthermore, is the number of discretized k points in the first Brillouin zone, and A^i 
and N2 determine the range of the Bloch bands to be included in the construction. In the 
case of NiO we select the five energy bands around the Fermi level, shown by the solid lines 
in Fig. [H for our present calculation, and the subscript n in enumerates the five Ni 3d 
LMTOs. These five Bloch bands are dominated by Ni 3d states, but also have considerable 
contributions from O 2p states resulting from the hybridization effect. The quantities |lV„fc) 
are simply the Fourier transformation of the real-space Wannier functions but in the 
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FIG. 1: Nonmagnetic band structure of NiO obtained by the LMTO method; the Fermi energy 
is set to zero. The five "d-hke" Bloch bands (in sohd hnes) are used for constructing Wannier 
functions. 



multi-band case they should be distinguished from eigenf unctions. Hereafter, we refer to 
\Wnk) as Wannier functions, too. 

So far we have a set of five nonorthogonal Wannier functions with d symmetry at every 
k point. To orthogonalize them, one needs to calculate the overlap matrix 0{k)nn' = 
{Wnk\Wn'k) and its inverse square root. The orthogonalized Wannier functions are given in 
the standard way, 

\Wnk) = J2\W^'k){0~Hk))^'n. (2) 
n' 

By using Eqs. ([T]) and and noticing that Hlipif.) = eikli^ik) and (0^|^/'ifc) = c„i(^) where 
ei{k) and c„j(fc) are the eigenvalues and eigenvectors of the single-particle Kohn-Sham equa- 
tions, one can express the LDA Hamiltonian matrix within the basis of orthogonalized Wan- 
nier functions as 

H'^^ikU' = {Wnk\H\Wr.,k) 

= Y.io{k)-'/\^ cak)emcHikmk)-'/'U'n'. (3) 

m,m' i=Ni 

From the above equation and the relation 0{k)nn' = J2i'=Ni ^ni{k)c*nn{k), it is easy to see 
that the matrix elements of H^^{k) can be built up exclusively in terms of the eigenvalues 
ei{k) and eigenvectors c„j(fc) of the LDA single-particle problem. Details of the procedure 



can be found in Ref. 52 



Next we present the density of states (DOS) of these Wannier functions in Fig. [2] and 
compare with the Ni 3d LMTO DOSs. The Wannier functions constructed in the above 
procedure have d symmetry, and therefore in a cubic system as NiO they can also be classified 
into threefold degenerate t2g states and twofold degenerate Cg states. However, one should 
keep in mind that these Wannier functions are not pure d states, but rather the anti-bonding- 
like states resulting from the hybridization between Ni 3d and O 2p states. This is obvious 
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FIG. 2: (color online). t2g and eg resolved LDA DOSs in the basis of d-like Wannier functions 
(t2g-red solid curve, eg-green dot-dashed curve) and Ni 3d LMTOs (t2g-black dotted dash curve, 
eg-blue dash curve). 



from Fig. [2], where in the Ni 3(i-dominant region the Wannier functions have more spectral 
weight than LMTOs, the extra spectral weight being due to the contribution of O 2p states. 
By construction, the bands at higher binding energy with mainly oxygen character are not 
included in the Wannier functions. Thus, as discussed above, the present scheme does not 
fully include the O 2p bands and the p-d hybridization. 



C. Inclusion of correlation effects by LDA-|-DMFT(QMC) 

The single-particle Hamiltonian matrix H^^{k) in Wannier function basis must be sup- 
plemented with the multi-orbital Coulomb interaction including the local Hund's rule ex- 
change interaction. In principle, i7^^(fc) should be corrected for the Coulomb interaction 
that has already been taken into account at the LDA level, yielding a double-counting-free 
single-particle Hamiltonian matrix H^^{k). But in the present case, this effect can be ap- 



proximated by a trivial shift of the energy (for details see Ref. l33l ). The full Hamiltonian 
thus reads 

k,n,n',cr i,n,n' ,cr,a' 

~ 2 •^nn'dlna^ln'a^in'a^ina- (4) 



The general interaction [/^^, denotes the Coulomb interaction matrix among different Wan- 
nier orbitals. For a cubic system it satisfies the relation = U — 2J(1 — 5„„/) — J5aa' to 
a good approximation. Here U corresponds to the effective interaction between electrons in 
Wannier orbitals as discussed above and not to the hare interaction between the actual d 
electrons. The values of the parameters U and J are accessible from a constrained LDA cal- 
culation. Specifically, from a constrained LDA calculation the average Coulomb interaction 
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U ^ and the Hund's rule exchange couphng J is obtained, which satisfies,— 

U + {M-l){U- 2J) + (M -!)([/ - 3J) 

^ = 2M^:i ' 

where M = 5 is the number of the interacting orbitals. Given U and J one can thus 
obtain U from Eq. ([5]). Since the constrained LDA calculation has not been performed for 
the Wannier functions, we take the U value for the LMTO basis from the literature,— and 
estimate the value for the Wannier function basis by 0"^^ = f/LMTO^]^ _ -^j^gj^g x is the 
admixture of oxygen 2p states into the Wannier functions.— For NiO, x is approximately 
0.15, and [/lmto _ g gy ii xhis means that f/^^ ^ 5.8 eV. Concerning the J value it is 
reasonable to assume that it is only weakly basis-dependent, and we just take the LMTO 
value J = 1 eV here. With these considerations we obtain U = 8 eV for the Wannier 
function basis. 

We solve the above material-specific model (jl]) by means of DMFT,— i^*^ which maps the 
original lattice model onto a single-impurity model subject to a self-consistent condition. The 
impurity problem is in turn solved by quantum Monte-Carlo (QMC) techniques,— yielding 
an imaginary-time Green function. The maximum-entropy method (MEM)^ is then used 
to obtain the physical spectral function on the real frequency axis. We neglect the spin-flip 
term in Eq. (j4]) which is known to lead to "minus-sign" problems in QMC simulations,— 
assuming that this term does not play a signiflcant role here. QMC simulations can only be 
performed at not too low temperatures, since the computation effort scales with 



III. RESULTS AND COMPARISON WITH EXPERIMENT 

Our DMFT(QMC) calculations were performed at temperatures T = 1160 K and 725 K 
respectively. In QMC simulations we used 40 time slices for the former case and 64 for the 
latter, and in both cases up to 10^ QMC sweeps were made to ensure good convergence. The 
electronic energy spectrum obtained at temperature T = 1160 K for interaction parameters 
U = 8 eV and J = 1 eV is shown in Fig. [31 The dominant peak at about —2.5 eV 
is due to the t2g bands which are hence completely fllled. The eg bands are split into 
lower and upper Hubbard bands. The insulating gap is therefore situated between the 
occupied t2g bands and unoccupied eg bands.— The shoulder between —6 eV and —4 eV 
originates from the occupied eg bands. We understand that this picture differs from the one 
provided by the cluster model,- where both the main peak and side peak are dominated 
by oxygen character. Indeed, in the present calculation, the hybridization between Ni 3d 
states and O 2p states is flxed at the LDA level, and a more complete theory of NiO should 
also allow the p-d hybridization to evolve under the influence of the interaction among d 
electrons. As discussed earlier this requires an explicit inclusion of oxygen states which is 
not implemented in the present work. Keeping the drawback of the present scheme in mind, 
we will leave the resolution of this problem to future studies, and take the present result 
as a flrst approximation. But as shown below, the agreement between the energy spectrum 
obtained in this calculation and the experiment is already very good. 

In Fig. m our LDA-I-DMFT results are compared with data of the combined XPS-BIS 
experiment by Sawatzky and Allen^ for the low-temperature AFM phase. Although our cal- 
culation was performed for the high-temperature PM phase a comparison is meaningful since 
the electronic structure of NiO is almost unaffected by the magnetic phase transition.— In 
view of the ambiguity in determining the position of the chemical potential in the insulating 
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FIG. 3: (color online). Theoretical energy spectrum of NiO obtained by the LDA+DMFT approach 
for T = 1160 K, and U = 8 eV, J = 1 eV. t2g states (blue dashed curve) and eg states (green 
dotted curve) are resolved, and total states (red solid curve) are the sum of the two. Note the 
curve for total states coincides with that of the unoccupied eg states above the chemical potential. 
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FIG. 4: (color online). Theoretical spectrum is obtained by the LDA+DMFT approach (red curve) 
for U = 8 eV and J = 1 eV at T = 1160 K, compared with experimental XPS+BIS data (black 
dots) after Sawatzky and Allen.— The zero energy point of the theoretical curve is shifted to fit the 
energy scale of the experimental data. 



gap we shifted the spectrum appropriately with respect to the experimental data. The over- 
all agreement between the calculated single-particle spectrum and the experimental data 
is surprisingly good. Only the shoulder below the main peak is not very well reproduced. 
It should be noted that the energy spectrum in this region, i.e., far away from the Fermi 
level, is difficult to obtain accurately due to the limitations imposed by MEM. On the other 
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FIG. 5: (color online). Theoretical spectrum of NiO obtained by the LDA+DMFT approach for 
U = 8 eV, J = 1 eV at T = 1160 K (red solid curve) and T = 725 K (green dashed curve) 
respectively. 

hand, matrix element effects, which have so far received little attention in the theoretical 
calculation of the XPS of NiO, may also play a role here. Although detailed quantitative 
information on matrix element effects is not available for NiO, it is generally known that in 
the XPS regime (i) O 2p features become weak relative to the Ni 3d features,— and (ii) the 
photoionization cross section of transition metals is suppressed at higher binding energies 
compared to lower ones.— The first effect is unlikely to be significant since oxygen fea- 
tures are only contained indirectly in the calculation. The second effect may even improve 
the agreement further by reducing the weight of the theoretical spectrum at higher binding 
energies. As a whole, we do not expect matrix element effects to impair the agreement of 
our results with experiment. 

To analyze the effect of the QMC simulation temperature on the spectrum we also per- 
formed calculations at T = 725 K. The result is shown in Fig. [5] in comparison with that 
obtained at T = 1160 K. Within the numerical error of QMC and MEM the spectra at 
those two temperatures do not show any significant difference. For this reason we expect 
our results to describe the spectrum of paramagnetic NiO also at even lower temperatures. 

It is also interesting to investigate how the spectrum changes when the interaction pa- 
rameter U is changed. For this purpose we performed the calculations also for U = 7.5 eV 
and 8.5 eV, with fixed J = 1 eV, at T = 1160 K. The results are shown in Fig. El in com- 
parison with that for U = 8 eV and the experimental data. Here, to illustrate the change of 
the insulating gap for different U values, the positions of the conduction-band peaks were 
kept fixed. From Fig. [6] we see that the best agreement with experiment is indeed found for 
U = 8 eV — the value determined from the constrained LDA calculation. The Mott gap is 
seen to increase roughly linearly with U. The t2g and eg splitting, however, does not change 
since it is controlled by the energy J whose value is fixed here. 

In Table |T] we present our results for the energy gap and the local magnetic moment in 
comparison with experimental results and the results of other theoretical approaches. Here 
we measured the energy gap of 4.3 eV as the distance between the half-maximum of the 
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FIG. 6: Theoretical spectra obtained by the LDA+DMFT approach for U = 7.5, 8, and 8.5 eV (blue 
dashed, red solid, and green dashed curves) respectively with fixed J = 1 eV and temperature T = 
1160 K. Experimental XPS+BIS data from Sawatzky and Allen^ are also shown for comparison. 

highest valence peak and the conduction peak, respectively, as was done in Ref. 0- The local 
magnetic moment of 1.70 fis is obtained in an indirect way. In experiment the local moment 
is associated with the spin polarization in the immediate vicinity of the Ni ion. However, 
in the DMFT calculation the magnetic moment is associated with the Wannier functions. 
For the parameter U = 8 eV and J = 1 eV, we find a magnetic moment associated with 
the Wannier functions as m «i 1.99 We can then estimate the local magnetic moment 

associated with the Ni ion by the relation mioc = 'fn{l — x) where x = 0.15 is again the 
contribution of O 2p states to Wannier functions.— This leads to mioc ~ 1-70 hb which is in 
accordance with the measurement and theoretical results obtained for the AFM phase (see 
Table [T]). Apparently the value of the local magnetic moment is almost unaffected by the 
AFM-PM phase transition. 



TABLE I: The LSDA, GW, LDA+U, LDA+DMFT (present work) and experimental energy gaps 
and magnetic moments for NiO. 





LSDA 


GW" 


LDA+U^ 


LDA+DMFT 


Expt. 


Energy gaps (eV) 


0.3 


3.7 


3.7 


4.3 


4.3^ 4.0^=^ 


Magnetic moments {^b) 


1.09 


1.83 


1.70 


1.70 


1.64^ 1.77^, 1.903 



""S. Massidda et alM. ^V. I. Anisimov et alM- 

""G. A. Sawatzky and J. W. Alleii^ '^S. Hiifner et al^ 

'^H. A. Alperin^ -^B. E. F. Fender et al^ 
9 A. K. Cheetham and D. A. O. Hope^ 
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IV. CONCLUSION 



In this paper we used LDA+DMFT to compute the electronic single-particle spectrum of 
NiO. Specifically, Wannier functions were constructed corresponding to the five Bloch bands 
across the Fermi level from the LDA band structure. These were used as the basis for a 
material-specific multi-band Hubbard-like Hamiltonian. This Hamiltonian was then solved 
by DMFT using QMC as the impurity solver. The electronic spectrum obtained in this way 
is in very good overall agreement with the experimental XPS+BIS results^ around the gap 
region. The energy gap and local magnetic moment obtained by us is in good agreement 
with experimental results. Taken together with the experimental observation that the long- 
range AFM order has no significant infiuence on the valence band photoemission spectra,— 
our result proves that the large insulating gap in NiO is due to strong electronic correlations 
in the paramagnetic phase. The magnetic order in NiO is therefore only a secondary effect, 
i.e., a consequence rather than the origin of the gap. Hence we expect the conduction band 
photoemission spectra to remain almost unchanged by the AFM long-range order. 

In the construction of the ab initio Hamiltonian (jl]), only the five "anti-bonding" bands 
(which have mainly Ni 3d characters in the LDA calculation) were included, while the three 
"bonding" bands (which are the mixture of the Ni 3d and O 2p states, but have more O 
contributions) below them were neglected. Therefore in the present work the contributions 
to the Wannier functions from Ni 3d and O 2p states depend on the LDA results. This 
means that the valence bands close to the Fermi level have mainly Ni 3d characters. The 
insulating gap we obtained is therefore an effective Mott-Hubbard gap between Wannier 
states. In other words, we used an effective Mott-Hubbard gap to mimic a charge-transfer 
gap^i^ or a gap with mixed character.— Strictly speaking, the gap here is also of mixed 
character since some amount of oxygen contribution is contained in the Wannier functions. 
The results obtained by such a treatment are surprisingly good. This may be due to the fact 
that correlation effects are treated better within DMFT than within any other theoretical 
approach available so far, and also because features due to oxygen are rather suppressed in 
XPS.-2S 

In spite of the limitations of the current implementation, we showed that the 
LDA-I-DMFT approach which combines first-principles, material-specific information with 
strong correlation effects is able to deal with late transition- metal monoxides like NiO. A 
more complete treatment of NiO within the LDA-I-DMFT approach will require the inclu- 
sion of oxygen bands and the p-d hybridization. Only in this way can one produce a full 
spectrum of NiO and identify the satellite structure at high binding energies. 
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